Engineered valley-orbit splittings in quantum confined nanostructures in silicon 
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An important challenge in silicon quantum electronics in the few electron regime is the poten- 
tially small energy gap between the ground and excited orbital states in 3D quantum confined 
nanostructures due to the multiple valley degeneracies of the conduction band present in silicon. 
Understanding the "valley-orbit" (VO) gap is essential for silicon qubits, as a large VO gap prevents 
leakage of the qubit states into a higher dimensional Hilbert space. The VO gap varies considerably 
depending on quantum confinement, and can be engineered by external electric fields. In this work 
we investigate VO splitting experimentally and theoretically in a range of confinement regimes. We 
report measurements of the VO splitting in silicon quantum dot and donor devices through excited 
state transport spectroscopy. These results are underpinned by large-scale atomistic tight-binding 
calculations involving over 1 million atoms to compute VO splittings as functions of electric fields, 
donor depths, and surface disorder. The results provide a comprehensive picture of the range of VO 
splittings that can be achieved through quantum engineering. 

PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.Gv, 71.70.Ej 



The ability to generate and manipulate three dimen- 
sionally confined quantum states in silicon down to the 
single electron regime is a much sought-after goal both 
in semiconductor quantum computing (QC) and quan- 
tum electronics. Silicon has not only been the primary 
platform of the semiconductor industry for over half a 
century, but in the quantum regime it also offers the ad- 
vantage of long spin coherence times [1] necessary for 
QC. Over a decade, steady progress has been made to- 
wards realizing quantum dot (QD) [H, Q and donor 
based single electron states to encode and process quan- 
tum information. However, non-trivial challenges to- 
wards quantum electronics in Si arise from the existence 
of multiple conduction band (CB) valley degeneracies in 
Si, which multiples the orbital degrees of freedom. Since 
qubits require a two-level spin system well isolated in en- 
ergy from other states, a critical goal is to establish con- 
trol over valley-orbit (VO) splitting by means of external 
perturbations and engineered quantum confinement. To 
date, a comprehensive understanding of VO splitting in 
Si has been hindered by a lack of consistent experimental 
data, with VO splittings measured over several orders of 
magnitude from ueV to meV 13l-[l7|. Although a num- 
ber of investigations made inroads towards a theoretical 
understanding of this problem [8-12|, a unified theory is 
still absent primarily due to the lack of realistic models of 
the region between Si and an insulator. In this Letter, we 
present measurements of the VO splittings in Si metal- 
oxide-semiconductor (MOS) nanostructures under vari- 



ous quantum confinement regimes engineered by internal 
and external electric fields. Million atom tight-binding 
(TB) simulations of the various confinement regimes are 
performed to explore the range of VO splittings possible 
in these structures, and to provide a unified theoretical 
underpinning of the experimental data. With this com- 
bined experimental and theoretical approach over a range 
of confinement regimes we are able to not only explain 
the range of VO splittings observed here, but in addition 
understand how to control the VO splitting in a range of 
device configurations for future applications. 

In general VO splittings are expected to depend criti- 
cally on the details of the confinement potential, interfa- 
cial disorder, barrier material, lattice miscuts, substrate 
orientation, strain, electric and magnetic fields ilii. 
Once understood, this suggests the ability to engineer 
VO interaction externally, and hence to directly tune the 
momentum space properties of Si. 

Hall-bar experiments in strained Si quantum well with 
SiGe barrier have reported VO splittings of the order 
of /ieVs |13l4l5|. It was shown that a strong vertical 



magnetic field, which reduces the lateral extent of the 
wavefunction and hence the exposure to disorder, could 
be used to obtain a relatively larger VO of 1.5 meV in 
the SiGe system [l5|. Relatively little data exits for VO 
splittings in Si MOS QDs near the few electron regime. 
In Ref [16], a VO splitting of 0.76 meV was reported 
recently. A past measurement of VO in a SIMOX device 
showed an unusually high value of 23 meV [l^, which is 
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yet to be explained conclusively, although a recent work 
has given some plausible arguments for the cause 0. 
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FIG. 1: a) SEM image of a FinFET device used in the exper- 
iments. b)-e) Schematic of the various confinement regimes 
showing the vertical confinement potential near the oxide in- 
terface, b) is a donor bound Coulomb confinement regime 
at low E- field, c) is a hybrid confinement regime between 
the donor and the interface well realized at higher E-fields 
when the two wells are lined up in energies and are strongly 
tunnel coupled, d) is an interfacial confinement regime real- 
ized at even higher fields, but laterally bound by the donor 
Coulomb potential, e) is a QD-like confinement regime re- 
alized at strong E-fields for device samples without any in- 
fluence of the donor. The lateral confinement is provided by 
the residual barriers in the access regions. The insets show 
schematics of the lateral confinement potential in d) and e). 

Fig. 1(a) shows an scanning electron micrograph 
(SEM) image of a MOS FinFET device used in the ex- 
periments reported here. A silicon nanowire connected 
to source and drain leads forms the channel of the tran- 
sistor. A second nanowire is deposited perpendicular to 
the channel as the gate electrode. A thin nitrided oxide 
layer separates the gate from the channel. 

This device has been used to realize and probe dif- 
ferent confinement regimes in Si nanostructures. The 
measurements are based on excited state transport spec- 
troscopy of a localized electron in a single quantum dot 
or donor. Most measurements utilize open system and 
investigate the quantum Hall effect in high B-field, and 
then extrapolate the measured value to the low B-field 
limit. Our method in essence provides an all-electrical 
means of measuring valley-orbit splitting. As noted, VO 
splittings change significantly with the applied electric 
field ranging from /ieVs to tens of me Vs. Hence, it is 
important to report the field at which VO measurements 
are done in experiments. Some device samples contained 
a single Arsenic donor in the channel due to diffusion 
from the leads, and could be used to probe donor states 
lol , [2q| . Other samples without donors could be used to 
realize QD type states. 

To provide a comprehensive theoretical basis for the 
experimental results we use a 10-band sp^d^s* nearest 



neighbor TB model involving over a million atoms, 
tailed in Refs 



de- 



23, M, The method has been used 



on Si nanostructures lo|, [H, [20|, [13, [29| . A dangling bond 
passivated surface model has been used to represent the 
interfacial boundary [30|, except when modeling the Si02 
insulator layer explicitly, as shown in Fig. 4. 

Fig. 1(b)- (e) show the schematic of the different con- 
finement regimes investigated in this work. The gate po- 
tential generates a triangular well near the oxide inter- 
face. If a donor is present in the channel, an additional 
Coulomb potential well forms on top of the triangular 
well at some distance from the interface, making it pos- 
sible to study different confinement regimes (Fig. 1 cap- 
tion) in a number of device samples. 

First, we discuss the three confinement regimes of Fig. 
1 (b), (c) and (d) that are influenced by the presence of 
a donor. In the experimental device, the As donors are 
located less than 6 nm from the oxide interface and are 
subjected to fields of tens of MV/m. This donor-interface 
configuration has been proposed as an important system 
for implementing a donor-dot hybrid qubit |23l -l' 
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FIG. 2: a) Tight-binding calculations of VO splitting in Si 
in the presence of a single As donor as a function of field 
for two donor depths. The inset shows a ID schematic of 
the confining potential with the donor, b) VO splitting as a 
function of donor depth and electric field obtained from TB. 
The white markers show the measured VO splitting for six 
device samples extracted from the measurements of Ref [l^ . 

For bulk donors at zero fields, all six valleys contribute 
to VO splittings. Since the central cell potential varies 
from one donor species to another, so does the VO split- 
ting [23I: in bulk the VO splitting of a bulk As donor is 
about 21 meV, compared to 12meV for a P donor. If an 
electric field is applied in the z direction, the weight of the 
wavefunction increases in the kz valleys and diminishes 
in the others. This wavefunction redistribution in mo- 
mentum space causes a reduced VO splitting. Fig. 2(a) 
shows the calculated VO splitting as a function of the 
E-field for two different donor depths. At low E-fields, 
the VO splitting is about 20 meV, comparable to the VO 
splitting of a bulk As donor. As the E-field is increased, 
the donor states hybridize with interface states, and VO 
splitting reduces gradually. At high enough E-fields, the 
electron is pulled to the interface, reducing VO splitting 
to a few meV, as expected of QD bound states at strong 
fields. Once the electron resides at the interface, the VO 
splitting varies linearly with the field, as shown by the 
red curve of a donor at 5.4 nm depth at fields above 30 
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MV/m. The blue curve is for a donor at a shallower 
depth of 2.7 nm. The change in VO is smoother because 
of stronger tunnel coupling between the two wells. 

In Fig. 2(b), the VO splitting is plotted in color code 
as a function of donor depth and electric field. At low 
fields, the donor bound Coulomb-like regime of Fig. 1(b) 
is realized, whereas at high fields, the states are mostly 
interface bound, as in Fig. 1(d). The measured VO split- 
tings for six devices with donors at various depths and 
fields are mapped on this figure as white markers. The 
data points sample out all three confinement regimes. 
The measured VO data has been extracted from mea- 
surements of DO excited state spectroscopy as reported 
in Ref [l9|, where the donor depths and applied fields 
were reported for various device samples. It is therefore 
possible to obtain a whole range of VO splitting in single 
donor devices ranging from 20 meV to an meV or even 
less. This means that VO splittings can be engineered 
through donor implantation depths and applied E-fields. 
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FIG. 3: a) Conductance vs. gate and bias voltage plot for a 
device without the influence of a donor (confinement regime 
Fig. 1(e)), showing blocked diamond region and tunneling 
through the QD states, characteristic of Coulomb blockade, 
b) Inset: Measured conductance data of the device as a func- 
tion of gate voltage and temperature (T), between 1.5 and 90 
K. The temperature steps are not all equal. Main plot: The 
peak maxima (extracted from the Inset) vs. 1/T, showing a 
cross-over from single level transport (linear regime) to classi- 
cal transport at higher T. c) Calculated spectrum of a MOS 
QD relative to the ground state, with lateral dimensions of 
30.5 nm x 38.5 nm, as a function of the vertical electric field, 
showing VO split states, matched to the measured values. 

In Fig. 3, new measurements are shown for a device 
with a QD-like confinement regime described in Fig. 1(e). 
Unlike the previous cases, the lateral confinement here is 
provided by the residual barriers in the access regions. 
Hence, these are more extended interface bound states, 
as expected of MOS QDs. The measured experimental 
data for transport through the gate-field confined states 



is shown in Fig. 3(a). The low charging energy of about 
10 meV, determined by the height of the Coulomb block- 
ade diamonds, indicates there are no dopants present. 
Typical charging energies for dopant bound states range 
from 30 to 50 meVs [l9[. Furthermore, it was found 
that spin filling is consistent with the first diamond cor- 
responding to the the first electron [21 1. 

To make sure there is no unobserved low lying state 
between the ground and the first excited state in the sta- 
bility diagram, we measured the temperature dependence 
of the low bias trace versus gate voltage. This method 
gives us the position of the lowest excited state [22| . The 
inset of Fig. 3(b) shows the temperature dependence of 
the traces. In the temperature regime where the peak 
conductance is a linear function of the inverse tempera- 
ture, 1/T, there is transport through a single quantum 
state. Here the peak conductance increases with decreas- 
ing temperature (Fig. 3(b)). At higher temperatures, 
multilevel based classical transport mechanisms domi- 
nate, and at even higher temperatures, thermally acti- 
vated transport dominates. The position of the cross-over 
point between the quantum regime and classical regime 
is determined by the position of the first excited state. 
Here, we found VO to be - 12 K, or - 1 ± 0.1 meV. 

A variety of factors such as orbital confinement, valley 
polarization, and VO interaction compete to determine 
the electronic structure of Si QDs. We have therefore 
performed TB simulations for this QD-like confinement 
regime as shown in Fig. 3(c). Since the effective lat- 
eral dimensions and the microscopic E-fields (F) in the 
channel are not known, we considered a range of val- 
ues for these parameters, and obtained best agreement 
with the measured data for lateral confinement lengths 
of 30.5 nm, Ly 38.5 nm and F 11.6 MV/m. 

An instance of a simulated energy spectrum is shown in 
Fig. 3(c) as a function of the E-field. The energy of the 
excited states are plotted relative to the ground state en- 
ergy. The orbital states of the dot arising from the same 
valley configuration appear as flat lines with the vertical 
E-field. Their energies are primarily determined by the 
lateral dimensions of the box. The shorter the confine- 
ment lengths, the higher are the orbital energies, con- 
sistent with a 2D particle in a box model. The VO split 
states appear as tilted lines, showing a linear dependence 
on the E-field. An 1 meV VO splitting is obtained at 
F = 11.6 MV/m independent of the lateral dimensions. 
The lowest three measured energy gaps are superimposed 
on this plot as white markers. To check that the 1 meV 
state is not an orbital state of the same valley configura- 
tion, we performed a back-of-the-envelope calculation to 
estimate the confinement length it would correspond to 
if it were indeed an orbital state, and obtained L = 77 
nm. Comparing with the channel and the gate lengths, 
the QD states are expected to have confinement lengths 
of less than 60 nm. Therefore, the 1 meV state is indeed 
a VO split state. 
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Since the VO split states can cross the orbital states 
of the dot at higher E-fields, the level arrangements in Si 
dots can vary from one experiment to another depending 
on the E- field and effective lateral confinement lengths. 
Hence, it is important to determine the E-field when re- 
porting a VO measurement. The first excited state en- 
ergy gap of the QD changes slope with E-field at this 
crossing point, which can be used as a reference point 
for measuring VO and orbital splittings in experiments 
in which the vertical field can be tuned. 
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FIG. 4: (c) TB simulations of VO splitting as a function of 
the vertical field for an ordered (a) and disordered (b) surface. 

VO splittings in realistic systems are likely to be influ- 
enced by the atomistic details of the Si02 — Si interface. 
We performed TB calculations to investigate the role of 
interface disorder on VO splittings in MOS QDs. We 
have used a virtual crystal (VC) 4 band sp^ TB model 
of Si02 [m . Since the oxide is amorphous or highly dis- 
ordered in reality, this is an approximation. However, 
given that ab-initio studies have shown some crystalline 
structure in the oxide near the interface [sH where the ex- 
ponential tail of the QD wavefunctions reside, this model 
is expected to provide a good qualitative picture. 

In Fig. 4(a), we included 2 nm width of a VC Si02 in 
the simulation domain along with Si. To simulate disor- 
der, we have used an additional unit cell of Si02 at the 
the Si02 — Si interface (Fig. 4(b)), and replaced some of 
the virtual Si02 atoms with Si randomly. This in effect 
creates a local variation in AEc, and an overall decrease 
in the average AEc at the interface. Simulations of 10 
randomly generated samples show a slightly decreased 
VO splitting as a function of field. The standard devi- 
ation of VO splittings is represented by the error bars. 
Thus interface disorder does not have a significant effect 
in the calculated VO splittings in this model. In the field 
regimes investigated, the vertical E-field creates the most 
significant change in VO splittings. 

We have investigated experimentally and theoretically 
the valley-orbit splitting in Si MOS devices under a range 
of quantum confinement conditions and provide a unified 
description of the large range of VO gaps observed in 
QDs and in hybrid single donor states. Large scale atom- 
istic tight-binding simulations confirm quantitatively the 
range of observed VO splittings, and shed light on the 
role of interfaces, disorder, and electric fields. In Si MOS 



QD states, VO splittings are shown to increase with the 
vertical E-field, and can be influenced by atomistic dis- 
order at low E-fields. Presence of single donors in these 
devices can yield VO splittings from 20 meVs to sub 
meVs depending on the field. This work enhances our 
understanding of VO interaction induced energy gaps in 
Si nanostructures, a critical consideration for Si qubits. 
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